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Abstract 

Background: Efflux is a widespread mechanism of reversible drug resistance in bacteria that can be triggered by 
environmental stressors, including many classes of drugs. While such chemicals when used alone are typically toxic 
to the cell, they can also induce the efflux of a broad range of agents and may therefore prove beneficial to cells in 
the presence of multiple stressors. The cellular response to a combination of such chemical stressors may be 
governed by a trade-off between the fitness costs due to drug toxicity and benefits mediated by inducible systems. 
Unfortunately, disentangling the cost-benefit interplay using measurements of bacterial growth in response to the 
competing effects of the drugs is not possible without the support of a theoretical framework. 

Results: Here, we use the well-studied multiple antibiotic resistance (MAR) system in E. coli to experimentally 
characterize the trade-off between drug toxicity ("cost") and drug-induced resistance ("benefit") mediated by efflux 
pumps. Specifically, we show that the combined effects of a MAR-inducing drug and an antibiotic are governed by 
a superposition of cost and benefit functions that govern these trade-offs. We find that this superposition holds for 
all drug concentrations, and it therefore allows us to describe the full dose-response diagram for a drug pair using 
simpler cost and benefit functions. Moreover, this framework predicts the existence of optimal growth at a non- 
trivial concentration of inducer. We demonstrate that optimal growth does not coincide with maximum induction 
of the mor promoter, but instead results from the interplay between drug toxicity and mar induction. Finally, we 
derived and experimentally validated a general phase diagram highlighting the role of these opposing effects in 
shaping the interaction between two drugs. 

Conclusions: Our analysis provides a quantitative description of the MAR system and highlights the trade-off 
between inducible resistance and the toxicity of the inducing agent in a multi-component environment. The results 
provide a predictive framework for the combined effects of drug toxicity and induction of the MAR system that are 
usually masked by bulk measurements of bacterial growth. The framework may also be useful for identifying 
optimal growth conditions in more general systems where combinations of environmental cues contribute to both 
transient resistance and toxicity. 



Background transient resistance does not rely on mutations, but can 

The resistance of bacteria to antibiotics has prompted in- be induced by a large class of chemicals commonly 

tense scientific research in the last several decades be- found in drugs (e.g. antibiotics and painkillers) and food 

cause it directly underlies the clinical treatment of preservatives. These chemicals are typically toxic to the 

infections [1]. While a large number of studies have fo- cell when used alone, but they can also induce resistance 

cused on mutation-driven resistance, recent attention to a broad range of agents. Consequently, they may 

has also shifted to transient, or "inducible", drug resist- prove beneficial to cells in the presence of multiple stres- 

ance taking place on much shorter timescales [2-9]. This sors. The net effect of a combination of chemical stres- 
sors can therefore be counterintuitive, because it is 
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and drug toxicity. Such situations may arise, for example, 
in the human digestive tract, where bacteria face a 
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cocktail of diverse chemical cues. The combined effects 
of multiple stressors, in general, have been studied for 
many decades in hopes of optimizing the clinical efficacy 
of combinatorial therapies [10-12]. More recently, the 
effects of drug interactions on the evolution of irrevers- 
ible (mutation-driven) drug resistance have also been 
recognized [13-16]. Drugs that interact synergistically to 
produce a strong toxic effect can accelerate the acquisition 
of mutations conferring drug resistance [14]. On the other 
hand, antagonistic drug pairs produce a weaker toxic effect 
but can slow the acquisition of resistance [13,15]. These 
results demonstrate an inherent trade-off between the tox- 
icity of the drug combination and its potential to facili- 
tate drug resistance [17]. They also raise an interesting 
question: do trade-offs between drug toxicity and resist- 
ance also play a role in inducible drug resistance? Fur- 
thermore, when cells are exposed a combination of toxic 
agents that potentially induce transient resistance, how 
are these trade-offs related to synergy or antagonism be- 
tween the given agents? 

To address these questions, here we study inducible re- 
sistance mediated by the MAR (multiple antibiotic re- 
sistance) system. The MAR system, present in many 
bacterial species, consists of an operon that confers 
efflux-mediated [18-21] resistance to a broad range of anti- 
biotics and can be activated by a host of chemical agents, 
including analgesics and food preservatives [2-5,7-9,22- 
25]. For example, in E. coli, high concentrations (5 mM) of 
salicylate, a well-known analgesic, lead to an increase in 
the minimum inhibitory concentration (MIC) of several 
antibiotics, and this effect has been partially attributed to 
the MAR system [4,26,27]. However, MIC alone masks 
the distinct effects of salicylate toxicity, antibiotic tox- 
icity, and induction of the MAR system (the "benefit" of 
salicylate) on the growth, because it provides only an ag- 
gregate measure of these effects. Therefore, the under- 
lying nature of the cost-benefit interplay remains hidden. 
Here we combine experiments with a simple phenom- 
enological model to fully characterize a range of effects 
resulting from the trade-off between drug toxicity and 
the benefit mediated by efflux pumps. 

Results 

Experimental characterization of salicylate-induced 
antibiotic resistance 

To study the trade-off between drug toxicity and 
induced resistance, we first measured the effects of sali- 
cylate and two protein synthesis inhibitors, chloram- 
phenicol and tetracycline, on cell growth. As expected, 
the effect of each drug alone is to slow cell growth as 
concentration increases. To quantitatively characterize 
this effect, we define growth cost as the reduction in 
growth rate of cells treated with one drug relative to the 
growth rate of untreated cells. In general, we find that 



cost functions are well-described by Hill functions with 
I<i, the concentration of drug i at which cost is half 
maximal, and n, the Hill coefficient (Figure la). This 
mathematical form is consistent with standard dose-re- 
sponse models [10]. In the presence of a high concen- 
tration of chloramphenicol or tetracycline, however, we 
find that adding salicylate can increase growth (Fig- 
ure lb, Additional file 1: Figure SI). This effect was 
previously reported for a single concentration of sali- 
cylate [26]; here we provide the entire non-monotonic 
dependence of this suppressive effect on concentration, 
allowing us to tease apart contributions of salicylates 
cost from its benefit (below). We also performed detailed 
quantitative measurements of another suppressive inter- 
action between an antibiotic, chloramphenicol, but with 
a different inducer, sodium benzoate (Additional file 1: 
Figure SI). 

Phenomenological cost-benefit model for MAR-induced 
drug resistance 

To quantitatively model the interplay between inducer 
cost and benefit in a two-drug environment that 
includes an inducer, we assume that the effects of the 
two drugs, in the absence of MAR induction, are inde- 
pendent in the sense that the relative growth rate of 
cells in the presence of both drugs (g SA ) equals the 
product of individual relative growth rates of cells in 
the presence of each drug alone (gsA = gs gA)- This as- 
sumption, known as Bliss independence [11], provides 
exact results for simple situations, such as a single en- 
zymatic reaction disrupted by mutually non- exclusive 
inhibitors [10], in which case the effects of the drugs 
are defined in terms of reaction fluxes. However, the 
assumption of Bliss independence often fails for more 
complex systems, and it is therefore primarily used as a 
phenomenological null model for the physiological 
effects of non-interacting drugs on cell proliferation or 
growth [10,14,28-30]. 

Here, we extend the concept of Bliss independence to in- 
clude drug interactions mediated by the induction of the 
MAR system. Specifically, we assume that the presence of 
one drug (S) provides a fitness benefit by reducing the ef- 
fective concentration of the second drug (A), thereby coup- 
ling the effects of the drugs. We incorporate this rescaling 
and re-write the model in terms of the combined growth 
cost, C S a (by definition 1- g SA ) to arrive at 

C SA =f(A eff ) +g(s) -f(A eff )g(s) (1) 

where f(A) and g(S) are the growth costs of drugs alone 
(Figure la). Equation 1 expresses Bliss independence in 
terms of growth costs, and also generalizes it to include 
an S -dependent reduction of A to A eff . Therefore, the 
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Figure 1 Single and Multiple Drug Growth Costs, a. Growth cost is defined as the reduction in growth rate of cells treated with a drug relative 
to the growth rate of untreated cells. Solid lines, best fits to Hill functions h(x) = x7(KP + x n ), with K, equal to the concentration of drug i (i = S for 
salicylate, i = A for both antibiotics) at which growth is inhibited by 50%. K A = 1.80 (1.66, 1.93) ug/mL, n = 1.97 (1.67, 2.26) for chloramphenicol; 
K A = 0.41 (0.37, 0.45) uoj/mL, n = 1 .90 (1 .61 , 2.1 9) for tetracycline; and K s = 6.06 (5.50, 6.62) mM, n = 1 for salicylate. 95% confidence intervals from 
nonlinear least squares fitting in parentheses. Error bars are standard deviation of four replicates, b. In the presence of high concentrations of 
chloramphenicol (Cm), growth is maximal for a nonzero concentration of salicylate (upper insets: OD time series). Contours of constant growth in 
concentration space indicate suppression. Dashed line, contours of 45% and 55% maximum growth. 
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Figure 2 Cost-Benefit Model for Interactions between MAR inducers and Antibiotics, a. The combined growth cost of two drugs (S and A) 
is assumed to be Bliss independent, as long as the effective concentration of one drug is reduced by the presence of the second drug. Thus, the 
combined cost of both drugs is approximately equal to the sum of effective growth costs, f and g, for each individual drug. The reduction of each 
drug's growth cost depends on the concentration of the other drug, which in turn dictates the growth benefit provided, b. mar promoter activity, 
defined as the steady state time-averaged Fluorescence/OD times relative growth rate, increases as a function of salicylate concentration. Error 
bars are +/- one estimated standard deviation. Inset: mar fluorescence concentration (YFP fluorescence/optical density) time series for 
concentrations of salicylate ranging from 0 mM (purple) to 6 mM (blue). Data points are means over two replicates. Promoter activity is estimated 
by averaging fluorescence concentration in the steady state and multiplying by relative growth rate. 

V ) 



approximate additivity of drug costs implied by Bliss in- 
dependence is modified to an approximate additivity of 
effective costs (Figure 2a). The model assumes that the 
presence of inducer reduces the effective concentration of 
drug A and adds an additional cost, but otherwise does 



not affect growth rate. This assumption considerably lim- 
its the spectrum of possible cellular responses to the drug 
pair, because drug S can only change the effective con- 
centration of drug A, but will not change the shape of its 
cost function. The model can be readily extended to 
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include two drugs which both act as inducers, in which 
case S is also changed to S eff . However, because the anti- 
biotics used are poor inducers of the MAR system 
(Additional file 1: Figure S7), we neglect their effects on 
the concentration S and limit this study to the direc- 
tional model implied by Equation [1]. 

A simple model, which assumes that the concentra- 
tion of antibiotic A is reduced by the induction of efflux 
pumps (Additional file 1), suggests a form for A eff : 



hypothesis, we can rearrange equation 1 to express this 
concept as 



(2) 



We call the function (3(S) the inducible benefit, and 
it contains all quantitative information about the re- 
sistance mechanisms induced by drug S. To estimate 
the inducible benefit, (3(S), for the MAR system, we 
measured the activity of the mar promoter as a func- 
tion of S (Figure 2b; Additional file 1). MAR promoter 
activity could dictate the functional form for (3(S) if, 
for example, we assume it is proportional to the 
number of efflux pumps produced in response to sa- 
licylate (Additional file 1). We scale this promoter 
activity by an adjustable parameter (3 max , so that 
y£( s ) =£ max S/(Ki n d + S), where K ind = 0.8 is the con- 
centration of inducer that yields half maximal promoter 
activity (Figure 2b). The adjustable parameter (3 max links 
mar activity induced by drug S to a phenotypic response 
(resistance) to a second drug A. This parameter will 
be specific to drug A, and it provides a measure of 
the efficiency with which MAR induction eliminates A 
from the cell 

The cost-benefit theory assumes that the effects of 
the inducing drug and the antibiotic are independent, 
up to a reduction in the concentration of the antibiotic. 
The aim of our model is not to achieve a microscopic 
theoretical description of the system, but rather to pro- 
vide a minimal phenomenological model that quantita- 
tively captures the measured behavior. The success or 
failure of the model must therefore be determined ex- 
perimentally. Specifically, our model does not attempt 
to elucidate the microscopic variables governing the 
multi-drug effects — which would require dozens, if not 
hundreds, of microscopic parameters — but rather posits 
that a simple relationship should exist between cell 
growth in the presence of one drug (which is a function 
of the cells internal state) and cell growth in the pres- 
ence of two drugs (which is a function of an entirely 
different intracellular state). Specifically, the model 
requires equality between cell growth in the presence 
of A and cell growth in the presence of the com- 
bination S and A (with A'>A), once we account 
for the costs of drug S. To directly verify this 



Q 



SA 



g(s) 



l-g(S) 



eff) 



(3) 



The left hand side represents the "adjusted" multi-drug 
cost, once the toxic effects of S (cost) have been 
removed. The model implies that this adjusted cost of 
the drug combination, as a function of A eff , is function- 
ally equivalent to the cost function f(x) for a single drug. 
Hence, if one removes from the multi-drug costs (C S a) 
the effects of g(S) according to the left-hand side of 
Equation 3 and then properly chooses the concentration 
reduction A^A eff (determined by the single parameter 
Pmax)> all two-drug data from a given drug combination 
should collapse to the single curve determined by f(x). 
Note that because a Hill function describes costs for 
many individual drugs as well as slices through many 
drug combination effect surfaces [10], collapsing the 
multi-drug cost curves from a given drug pair (Figure 3a) 
to a common Hill form can be simply achieved using in- 
dependent parameters for each salicylate concentration. 
By contrast, equation 3 suggests that a non-trivial col- 
lapse of all curves for a given drug pair requires only a 
single parameter (3 max . Once this parameter is deter- 
mined, the beneficial effects of different concentrations 
of salicylate are linked by the mar promoter activity 
curve. We verified this constraint (Figure 3b) by fitting 
the single free parameter, (3 max , using two-drug cost curves 
for salicylate and chloramphenicol (p max =1.15 +/- 0.15) 
and salicylate and tetracycline (p max = 6.04 +/- 0.16). In all 
cases, the model provides an excellent fit to the data 
(R 2 = 0.98; see Additional file 1: Figure S2 for a direct com- 
parison of growth rates from experiment and model). 
Interestingly, this analysis demonstrates that a superpos- 
ition of cost and benefit functions quantitatively describes 
the combined effects of salicylate and an antibiotic. Sur- 
prisingly, we find that the benefit function depends only 
on the inducer concentration and is dictated by the mar 
induction curve. 

Transition from salicylate as toxic to salicylate as 
beneficial 

Our results uncover an interesting range of cellular be- 
havior that emerges from the trade-offs of salicylate 
toxicity and the simultaneous induction of multi-drug 
resistance. First, for each concentration of antibiotic, we 
find that growth is maximal at a single salicylate con- 
centration S*. In addition, we see an apparent transition 
between two regimes — one where the presence of the 
inducing drug is harmful, and another where it is bene- 
ficial—as A eclipses a threshold A crit (Figure 4a). That 
is, S* becomes non-zero as A crosses A crit . In other 
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Figure 3 Cost-Benefit Analysis Quantitatively Describes Multi-Drug Growth Cost Functions, a. Multi-drug growth cost curves for salicylate 

and tetracycline (top) and salicylate and chloramphenicol (bottom) are increasing functions of antibiotic concentration, but also depend on the 

concentration of salicylate. Each growth cost curve represents growth in increasing concentrations of antibiotic at a single salicylate 

concentration, ranging from 0 to 7 mM. Black lines, ([salicylate] = 0 mM), which correspond to the single drug growth cost functions. Data points, 

means of four replicates. Error bars, +/- one sample standard deviation, b. The cost-benefit model requires rescaled versions of the multi-drug 

cost functions to have the same mathematical form as single drug cost functions (Equation 3). To directly test this assumption, we rescaled the 

multi-drug cost functions in panel a by removing the contribution from salicylate toxicity (cost) to each multi-drug growth cost function. Then, 

we rescaled the concentration of each antibiotic (chloramphenicol or tetracycline) using the single parameter (3 max , which accounts for the 

salicylate-induced reduction of intracellular antibiotic concentration (benefit). We trivially exploit the common Hill form of the single drug costs 

functions to show data from both drug pairs on a single plot, achieved by replacing drug concentration with effective toxicity, defined as ([A]/K A ) n , 

where K A and n characterize the single drug growth cost function for the antibiotics (see Figure 1). We use the salicylate induction curve of the mar 

promoter (Figure 2b) as the benefit function (3(S), for both drug interactions, with only the overall scale of the benefit function ((3 max ) specific to the 

two antibiotics. See also Additional file 1: Figure S2 for a direct comparison of growth rates from experiments and models. 
^ J 



words, the cell benefits from the presence of salicylate, 
but only when the antibiotic environment is sufficiently 
toxic to offset the inherent toxicity of the inducer. 
Interestingly, we find that mar promoter activity satu- 
rates at S^4 mM (Figure 2b), which is greater than S* 
for all concentrations of antibiotics used. Therefore, 
maximum suppression at S* results from a non-trivial 
interplay between MAR expression and inducer cost, 
not merely a maximum in MAR expression. 

Contributions of inducer cost and benefit to apparent 
MIC's of antibiotics 

Given our characterization of salicylate-induced multi- 
drug resistance, it is straightforward to calculate the ap- 
parent MIC of an antibiotic in the presence of any con- 
centration of salicylate. Here, we define the apparent 
MIC to be the minimum concentration of drug A at 
which relative growth is reduced to a value of 5. For 
small concentrations of salicylate, the apparent MIC's of 



both tetracycline and chloramphenicol increase dramat- 
ically (Figure 4b). However, at higher concentrations, 
the toxicity of salicylate begins to overwhelm its poten- 
tial benefits, eventually leading to a decrease in appar- 
ent MIC. It is therefore clear that MIC reflects 
contributions of both the cost and benefits of the indu- 
cing drug. By contrast, the interdependence of inducible 
resistance and mar promoter activity can be described 
using only a single parameter, (3 max , which provides a 
quantitative measure of inducible benefit that is not 
masked by the effects of inducer cost. 

Phase diagram for MAR-mediated drug interactions 

The interactions between salicylate and tetracycline, sali- 
cylate and chloramphenicol, and sodium benzoate and 
tetracycline are all suppressive, but this class of models 
describes a range of interactions ranging from synergistic 
to suppressive, based on the interplay of induction bene- 
fit and drug toxicity. Using Equations 1 and 2, we 
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Figure 4 Phase Diagram for Salicylate-Chloramphenicol Interactions, a. The concentration of salicylate, S* (in mM), at which growth is 
maximal depends on the external concentration of antibiotic. A transition occurs between the low antibiotic regime, where the presence of 
salicylate is harmful, to the high antibiotic where the presence of salicylate increases the growth of cells. Red, tetracycline; Blue, chloramphenicol. 
Concentrations for both drugs measured in ug/mL. Curves, numerical calculation from model; circles, estimates of S* from experiments at single 
concentrations of A. Maxima were estimated using cubic spline interpolation between data points at different concentrations of salcilyate. Error 
bars, concentrations corresponding to +/- 1% of maximum, b. Apparent MIC, defined as the minimum concentration of antibiotic (red, 
tetracycline; blue, chloramphenicol; both in units of ug/mL) at which growth is reduced to 6 = 0.5. While the quantitative results depend on the 
precise definition of MIC, the qualitative features do not depend on the choice of 6. c. Cost-benefit theory predicts general properties of drug 
interactions as a function of the maximum inducible benefit (3 max and the concentration of inducer, S. Solid line, phase boundary between 
synergistic and antagonistic drug interactions. Dashed line, phase boundary between antagonistic and suppressive interactions. In general, a 
higher value of the (3 max increases the antagonism between drug pairs at a given concentration of drug S. Parameters characterizing individual 
drugs include K s and K ind , which characterize, respectively, the cost of drug S and the corresponding induction of resistance systems (e.g. efflux 
pumps), and K A and n, which characterize the cost of drug A. By contrast, (3 max couples the individual effects of two drugs. Insets, heat maps of 
two-dimensional growth surfaces for 3 cell strains in the presence of Salicylate (Sal) and Chloramphenicol (Cm); top: (3 max = 1.15 (WT cells; 
suppressive), (3 max = 0.19 (mar mutant, antagonistic), and (3 max = —0.1 5 (tolC mutant, synergistic). See also Additional file 1: Figure S3. 
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quantitatively determine a phase diagram (Figure 4, Add- 
itional file 1: Figure S3) that specifies how a drug inter- 
action depends on the trade-offs between inducible 
benefit (described by the parameters K inc i and |3 max ), the 
cost of inducer (K s ), and the steepness of the antibiotic 
dose-response curve (n). Specifically, the cellular re- 
sponse to the first drug (S) alone includes both the 
growth cost of the drug (characterized by K s ) and the in- 
duction of MAR system (characterized by K in d). The dy- 
namics of the efflux pumps and their specificity for drug 
A determine the concentration reduction (A to A eff ), 
which is governed by p max . The phase diagram demon- 
strates that properties of the drugs alone (n, K s , Kind) de- 
termine the level of drug coupling, contained in (3 max , 
required to achieve antagonism or suppression. Drugs 
that strongly induce growth benefit and have low asso- 
ciated cost (K ind <<K s ) are always suppressive. By con- 
trast, high cost inducers (K inc i >>I<s) can never be 
suppressive because the cost of induction is too high 
(Figure 4c). 

The inducible benefit parameter (|3 max =1.15) charac- 
terizing the salicylate and chloramphenicol combination 
in wild type cells is far above the suppressive-antagonis- 
tic boundary ((3 max > > 2K ind /(K 1 n) = 0.45), and the inter- 
action is clearly suppressive for all concentrations of 
salicylate. To explore different regions of the phase dia- 
gram experimentally, we measured the effects of two dif- 
ferent mutations on the suppressive drug interaction 
between salicylate and chloramphenicol. The first mu- 
tant, AtolC [31], lacks the protein TolC required for ef- 
flux pumping [32]. While we expect that the AtolC 
deletion will partially suppress the benefit induced by sa- 
licylate, we cannot predict how the costs of the drugs 
and, therefore, the precise nature of the drug interaction, 
will be altered. Experimentally, this mutant showed 
decreased resistance to chloramphenicol (Figure 4a, 4b), 
corresponding to a rescaling of the single drug cost (i.e. 
K A smaller than in wild type), but the cost of salicylate 
remains unchanged (K s approximately same as wild 
type). |3 max was measured to be -0.15 +/- 0.01. Accord- 
ing to the phase diagram (Figure 4), the interaction 
should therefore be weakly synergistic, and in fact the 
contours of constant growth appear slightly concave 
(lower inset). The suppressive drug interaction has been 
eliminated because the benefit associated with salicylate 
has been decreased while its cost remains unchanged, 
resulting in significantly less antagonistic interaction. 

Since the suppressive interaction between salicylate 
and chloramphenicol is partially associated with the syn- 
thesis of AcrAB-TolC efflux pumps, we also hypothe- 
sized that cells with constitutive mar promoter activity 
would disrupt the cost-benefit interplay required for 
drug suppression. We selected such a mutant, here- 
called a tet mutant, by growing cells in 1 ug/mL 



tetracycline for 48 hours. We observed that the MAR 
system is a common target for mutations that confer re- 
sistance to tetracycline (see also [23]) (Additional file 1: 
Figure S6). In terms of our cost-benefit analysis, these cells 
synthesize resistance systems (benefit) without the asso- 
ciated cost of an inducing drug. To verify the molecular 
basis of resistance in the tet mutant, we sequenced the full 
genome and found only a single point mutation in the a5 
region of the MAR repressor marR, which is linked to 
dimer formation and subsequent binding to the mar pro- 
moter [33] . In the salicylate- chloramphenicol combination, 
the tet mutant showed increased resistance to chloram- 
phenicol (K A is larger than in wild-type) and a near- addi- 
tive drug interaction (p max = 0.19 +/- 0.03) between 
salicylate and chloramphenicol (Figure 4). Because the cost 
of salicylate is unchanged, the elimination of suppression 
suggests that the mutation has blocked the inducible bene- 
fit of salicylate (see also Additional file 1: Figure S5 for 
similar results with a different drug pair). Biologically, the 
inactivation of the MarR repressor no longer requires the 
presence of salicylate, and therefore the cost of the drug 
exceeds its associated inducible benefit. Interestingly, we 
also found that the mar activity in the mutant is higher in 
the absence of salicylate than the maximum mar activity 
induced by even high concentrations of salicylate in the 
wild type strain. However, adding salicylate further 
increases mar activity in the mutant (Additional file 1: 
Figure S7A). Moreover, the functional dependence of 
the induction on salicylate is similar to that of the wild- 
type (Additional file 1: Figure S7A, inset). This result is 
not intuitive because the beneficial effects of salicylate 
on the growth of the mutant in the presence of chloram- 
phenicol are quite small (p max is approximately 17% of 
that of the wild type, and the antagonism between the 
drugs is markedly decreased). 

Discussion 

Efflux is a widespread cellular defense mechanism with 
clinical consequences for antimicrobial treatments. Cost- 
benefit principles provide a convenient tool for analyzing 
these ubiquitous systems. We have shown that such a 
framework provides a quantitative description of drug 
interactions associated with the classic MAR system, 
which mediates an inducible resistance to many different 
toxic agents. The model allows us to quantify the relative 
importance of the induction mediating reversible resist- 
ance and the toxicity of the inducing agent. In particular, 
we observe a transition between two regimes, one where 
salicylate is harmful and one where it is beneficial, as the 
concentration of antibiotic crosses a threshold A crit . By 
disentangling the contributions of inducer cost and 
benefit, the analysis also provides a quantitative measure 
(Pmax) of the reversible (MAR-mediated) resistance to 
each antibiotic used. This measure complements 
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previous MIC measurements [26], which provide infor- 
mation about the net effect of both inducer costs and 
benefits but that mask the relative contributions of each. 
We demonstrate that it is precisely this balance between 
cost and benefit that determines whether the antimicro- 
bial properties of multiple agents are synergistic, antag- 
onistic, or suppressive. Therefore, in addition to 
providing a new quantitative model of the MAR system, 
we also lay the groundwork for a link between the short- 
term trade-offs of inducible resistance, which determine 
the type of interaction between pairs of toxic agents, and 
the acquisition of irreversible resistance mediated by 
these drug interactions. 

This framework may be applicable to other systems 
where environmental signals can be toxic but also induce 
transient resistance. For example, cancer cells can acti- 
vate pro-survival stress responses upon treatment with 
chemotherapies or radiation [34], and subpopulations of 
tumor cells can access new transcriptional/differenti- 
ation states to become reversibly drug resistant[35]. Our 
analysis implements efflux-mediated resistance using a 
concentration reduction A —> A eff for the effective con- 
centration of drug A, and this form could be used to de- 
scribe other common mechanisms of resistance, 
including enzymatic decay and target modification. As a 
consequence, the cost-benefit framework developed here 
may prove useful for modeling and understanding a wide 
range of systems governed by the inherent trade-offs of 
inducible resistance. 

Conclusions 

Our analysis provides a quantitative description of the 
MAR system and demonstrates that optimal cell growth 
in the presence of an inducer and an antibiotic does not 
coincide with maximum induction of the mar promoter, 
but instead results from the interplay between drug tox- 
icity and mar induction. In addition, we show that this 
interplay determines whether the drugs are synergistic, 
antagonistic, or suppressive. The framework may also be 
extensible to more general systems where combinations 
of environmental cues contribute to both transient resist- 
ance and toxicity. 

Methods 

Strains 

We used the E. coli Frag- IB strain, which is wild-type for 
the AcrAB multiple drug efflux pump system, and also 
the wild-type Exoli BW25113 strain. The AtolC mutant 
is strain JW5503-1 from the Keio collection. 

Drugs 

Drug solutions were made from solid stocks (sodium sa- 
licylate, Sigma-Aldrich; chloramphenicol, MP Biomedi- 
cals; tetracycline hydrochloride, Acros Organics; 



doxycycline hyclate, Sigma- Alder ich; kanamycin, Fisher; 
ciprofloxacin, Sigma- Alderich). All antibiotic stock solu- 
tions were stored in the dark at -20°C. All drugs were 
thawed (when necessary) and diluted in sterilized Lennox 
LB broth (Fisher) for experimental use. 

Growth conditions and drug treatments 

We inoculated LB media with a single colony and grew 
the cells overnight (12 h at 30°C with shaking at 
200 rpm). Following overnight growth, stationary phase 
cells were diluted (600-1000 fold) in LB media contain- 
ing 0.2% glucose and grown for an additional 3 h (30°C 
with shaking at 200 rpm). Following the initial dilution 
period, we transferred cells to 96 well microplates (round 
bottom, polystyrene, Corning) and set up a two-dimen- 
sional matrix of drug concentrations in each 96-well 
microplate (165 \A media per well). For the remainder of 
the experiment (4-8 h), cells were grown at 30°C with 
shaking at 1200 rpm on a Heidolph Titramax vibrating 
microplate shaker (Brinkmann). Plates were wrapped in 
aluminum foil to minimize evaporation. A600 (absorbance 
at 600 nm, proportional to optical density OD) and YFP 
fluorescence were measured at 15-25 min intervals for 
4-8 h using a Wallac Victor-2 1420 Multilabel Counter 
(PerkinElmer). 

Growth rate estimation and MIC determination 

Growth rates were determined by fitting multiple regions 
of growth curves (A 600 vs. time, Additional file 1: Figure 
S4) in early exponential phase to an exponential function 
(MATLAB 7.6.0 curve fitting toolbox, The Mathworks) 
using a variable length (45-150 minutes) sliding window. 
Growth rate corresponded to the window position asso- 
ciated with best fit to an exponential function, determined 
by a minimum root mean squared error (RMSE). For each 
experimental condition, the final growth rate was the sam- 
ple mean of 4 replicates (See also Additional file 1). 
Growth rate contours were estimated from experimental 
data using the Matlab csaps function to interpolate be- 
tween data points (typically 96) spaced roughly evenly 
throughout the two dimensional space of drug con- 
centrations. MIC was taken to be the drug concentra- 
tion^) at which growth was inhibited by 50%. 

Drug resistant mutants 

We diluted liquid cultures of wild-type cells (in station- 
ary phase) 1000-fold into 96 individual 150 \iL cultures 
on a single microplate, with each well supplemented with 
1.0 (ig/mL tetracycline. After approximately 48 h of 
growth at 30°C with shaking, optical density and fluores- 
cence were measured from each well. Samples from ran- 
domly selected cultures exhibiting high OD and 
fluorescence (Additional file 1: Figure S6) were frozen in 
15% glycerol at -80°C. Frozen cultures were used to 
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streak LB plates, and all experiments (including sequen- 
cing) were performed on cultures grown from a single 
colony isolated from these plates. 

Additional file 



Additional file 1: Supplemental Figures, Supplemental Methods, 
and Supplemental Notes [10,23,36-42]. 



Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

KW and PC designed research, analyzed data, and wrote the paper. KW 
performed all experiments. Both authors read and approved the final 
manuscript. 

Acknowledgements 

We thank L. Bruneaux, E. Balleza and A. Subramaniam for technical guidance, 
C. Guet for advice and for providing the Frag-1 B pZS*2 MAR-venus strain, 
and all members of the Cluzel lab for many helpful discussions. We also 
thank K. Dave for editorial advice and assistance and Kris Wood for 
comments on the manuscript. This work was supported in part by the NIH 
P50 award P50GM081 892-02 to the University of Chicago and the NSF 
Postdoctoral Fellowship 0805462 (to K.W.). 

Received: 22 December 201 1 Accepted: 24 April 2012 
Published: 25 May 2012 

References 

1. Levy SB, Marshall B: Antibacterial resistance worldwide: causes, challenges 
and responses. Nat Med 2004, 10(12 Suppl):S1 22-S1 29. 

2. Alekshun MN, Levy SB: The mar regulon: multiple resistance to antibiotics 
and other toxic chemicals. Trends Microbiol 1999, 7(10):41 0-41 3. 

3. Cohen SP, Hachler H, Levy SB: Genetic and functional analysis of the 
multiple antibiotic resistance (mar) locus in Escherichia coli. J Bacteriol 
1993, 175(5):1 484-1 492. 

4. Cohen SP, Levy SB, Foulds J, Rosner JL: Salicylate induction of antibiotic 
resistance in Escherichia coli: activation of the mar operon and a mar- 
independent pathway. J Bacteriol 1993, 175(24):7856-7862. 

5. Critzer FJ, Dsouza DH, Golden DA: Transcription analysis of stxl, marA, 
and eaeA genes in Escherichia coli 0157:1-17 treated with sodium 
benzoate. J Food Prot 2008, 71 (7):1 469-1 474. 

6. Herbert D, Paramasivan CN, Venkatesan P, Kubendiran G, Prabhakar R, 
Mitchison DA: Bactericidal action of ofloxacin, sulbactam-ampicillin, 
rifampin, and isoniazid on logarithmic- and stationary-phase cultures of 
Mycobacterium tuberculosis. Antimicrob Agents Chemother 1996, 40 
(10):2296-2299. 

7. Potenski CJ, Gandhi M, Matthews KR: Exposure of Salmonella Enteritidis to 
chlorine or food preservatives decreases [corrected] susceptibility to 
antibiotics. FEMS Microbiol Lett 2003, 220(2):181-186. 

8. Tavio MM, Vila J, Peril li M, Casanas LT, Macia L, Amicosante G, de Jimenez 
Anta MT: Enhanced active efflux, repression of porin synthesis and 
development of Mar phenotype by diazepam in two enterobacteria 
strains. J Med Microbiol 2004, 53(Pt 1 1):1 1 1 9-1 1 22. 

9. Viveiros M, Jesus A, Brito M, Leandro C, Martins M, Ordway D, Molnar AM, 
Molnar J, Amaral L: Inducement and reversal of tetracycline resistance in 
Escherichia coli K-12 and expression of proton gradient-dependent 
multidrug efflux pump genes. Antimicrob Agents Chemother 2005, 49 
(8):3578-3582. 

10. Greco WR, Bravo G, Parsons JC: The search for synergy: a critical review 
from a response surface perspective. Pharmacol Rev 1995, 47(2):33 1-385. 

1 1 . Bliss CI: The calculation of microbial assays. Bacteriol Rev 1 956, 20(4)243-258. 

12. Loewe S: The problem of synergism and antagonism of combined drugs. 

Arzneimittelforschung 1 953, 3(6):285-290. 

13. Chait R, Craney A, Kishony R: Antibiotic interactions that select against 
resistance. Nature 2007, 446(71 36):668-671. 



14. Hegreness M, Shoresh N, Damian D, Hartl D, Kishony R: Accelerated 
evolution of resistance in multidrug environments. Proc Natl Acad Sci USA 
2008, 105(37):1 3977-1 3981. 

15. Michel JB, Yeh PJ, Chait R, Moellering RC Jr, Kishony R: Drug interactions 
modulate the potential for evolution of resistance. Proc Natl Acad Sci USA 
2008, 105(39):1 491 8-1 4923. 

16. Chait R, Shrestha S, Shah AK, Michel JB, Kishony R: A differential drug 
screen for compounds that select against antibiotic resistance. PLoS One 
2010, 5(12):e15179. 

17. Torella JP, Chait R, Kishony R: Optimal drug synergy in antimicrobial 
treatments. PLoS Com put Biol 2010, 6(6):e 1000796. 

18. Levy SB: Active efflux mechanisms for antimicrobial resistance. Antimicrob 
Agents Chemother 1 992, 36(4):695-703. 

19. Levy SB: Active efflux, a common mechanism for biocide and antibiotic 
resistance. J Appl Microbiol 2002, 92(Suppl):65S-71 S. 

20. Okusu H, Ma D, Nikaido H: AcrAB efflux pump plays a major role in the 
antibiotic resistance phenotype of Escherichia coli multiple-antibiotic- 
resistance (Mar) mutants. J Bacteriol 1996, 1 78(1 ):306-308. 

21 . Wall ME, Markowitz DA, Rosner JL, Martin RG: Model of transcriptional 
activation by MarA in Escherichia coli. PLoS Comput Biol 2009, 5(12): 
e 10006 14. 

22. Cohen SP, Hooper DC, Wolfson JS, Souza KS, McMurry LM, Levy SB: 
Endogenous active efflux of norfloxacin in susceptible Escherichia coli. 
Antimicrob Agents Chemother 1 988, 32(8):1 1 87-1 1 91 . 

23. Cohen SP, McMurry LM, Hooper DC, Wolfson JS, Levy SB: Cross-resistance 
to fluoroquinolones in multiple-antibiotic-resistant (Mar) Escherichia coli 
selected by tetracycline or chloramphenicol: decreased drug 
accumulation associated with membrane changes in addition to OmpF 
reduction. Antimicrob Agents Chemother 1989, 33(8):1 318-1 325. 

24. Cohen SP, McMurry LM, Levy SB: marA locus causes decreased expression 
of OmpF porin in multiple-antibiotic-resistant (Mar) mutants of 
Escherichia coli. J Bacteriol 1988, 1 70(1 2):541 6-5422. 

25. McMurry LM, George AM, Levy SB: Active efflux of chloramphenicol in 
susceptible Escherichia coli strains and in multiple-antibiotic-resistant 
(Mar) mutants. Antimicrob Agents Chemother 1994, 38(3):542-546. 

26. Rosner JL: Nonheritable resistance to chloramphenicol and other 
antibiotics induced by salicylates and other chemotactic repellents in 
Escherichia coli K-12. Proc Natl Acad Sci USA 1985, 82(24):8771-8774. 

27. Price CT, Lee IR, Gustafson JE: The effects of salicylate on bacteria. Int J 
Biochem Cell Biol 2000, 32(1 0):1 029-1 043. 

28. Fitzgerald JB, Schoeberl B, Nielsen UB, Sorger PK: Systems biology and 
combination therapy in the quest for clinical efficacy. Nat Chem Biol 2006, 
2(9):458-466. 

29. Yeh P, Tschumi Al, Kishony R: Functional classification of drugs by 
properties of their pairwise interactions. Nat Genet 2006, 38(4)489-494. 

30. Lehar J, Zimmermann GR, Krueger AS, Molnar RA, Ledell JT, Heilbut AM, 
Short GF 3rd, Giusti LC, Nolan GP, Magid OA, et al: Chemical combination 
effects predict connectivity in biological systems. Mol Syst Biol 2007, 3:80. 

31 . Baba T, Ara T, Hasegawa M, Takai Y, Okumura Y, Baba M, Datsenko KA, 
Tomita M, Wanner BL, Mori H: Construction of Escherichia coli K-12 in- 
frame, single-gene knockout mutants: the Keio collection. Mol Syst Biol 
2006, 2:2006 0008. 

32. Fralick JA: Evidence that TolC is required for functioning of the Mar/ 
AcrAB efflux pump of Escherichia coli. J Bacteriol 1996, 1 78(1 9):5803-5805. 

33. Alekshun MN, Levy SB, Mealy TR, Seaton BA, Head JF: The crystal structure 
of MarR, a regulator of multiple antibiotic resistance, at 2.3 A resolution. 
Nat Struct Biol 2001, 8(8):71 0-71 4. 

34. Li F, Sethi G: Targeting transcription factor NF-kappaB to overcome 
chemoresistance and radioresistance in cancer therapy. Biochim Biophys 
Acta 2010, 1805(2):167-180. 

35. Sharma SV, Lee DY, Li B, Quinlan MP, Takahashi F, Maheswaran S, 
McDermott U, Azizian N, Zou L, Fischbach MA, et al: A chromatin-mediated 
reversible drug-tolerant state in cancer cell subpopulations. Cell 2010, 
141(1):69-80. 

36. George AM, Levy SB: Amplifiable resistance to tetracycline, 
chloramphenicol, and other antibiotics in Escherichia coli: involvement 
of a non-plasmid-determined efflux of tetracycline. J Bacteriol 1983, 

155:531-540. 

37. Alekshun MN, Levy SB: Regulation of chromosomally mediated multiple 
antibiotic resistance: the mar regulon. Antimicrob Agents Chemother 1997, 
41:2067-2075. 



Wood and Cluzel BMC Systems Biology 2012, 6:48 
http://www.biomedcentral.eom/1752-0509/6/48 



Page 11 of 1 1 



38. Maneewannakul K, Levy S: Identification of mar mutants among 
quinolone-resistant clinical isolates of Escherichia coli. Antimicrob Agents 
Chemother 1996, 40:1695-1698. 

39. Guet C, Bruneaux L, Oikonomou P, Cluzel P: In preparation.: ; 201 1:18. 

40. Le TT, Emonet T, Harlepp S, Guet CC, Cluzel P: Dynamical determinants of 
druginducible gene expression in a single bacterium. Biophys J 2006, 
90:3315-3321. 

41. Lutz R, Bujard H: Independent and tight regulation of transcriptional units 
in Escherichia coli via the LacR/O, the TetR/O and AraC/11 -12 regulatory 
elements. Nucleic Acids Res 1997, 25:1203-1210. 

42. Elf J, Nilsson K, Tenson T, Ehrenberg M: Bistable bacterial growth rate in 
response to antibiotics with low membrane permeability. Phys Rev Lett 
2006, 97:258104. 



doi:1 0.1 186/1 752-0509-6-48 

Cite this article as: Wood and Cluzel: Trade-offs between drug toxicity 
and benefit in the multi-antibiotic resistance system underlie optimal 
growth of E. coli. BMC Systems Biology 201 2 6:48. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at r~\ oz^wLm^ rpntral 

www.biomedcentral.com/submit \^ ™««™" ^enirai 



